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' Abstract 

^ ■ We show that superheavy threshold corrections in the New Minimal Supersym- 

metric GUT based on the SO(IO) Higgs system 210 126 126 10 120 can 
comfortably correct the prediction for the value of az{Mz) from the relatively large 
^ ' value predicted by the two loop RG equations to the central value determined by the 

0^ . current world average. The unification scale is raised above the one loop value over 

^ \ almost all of the viable parameter space. 
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1 Introduction 



Since the discovery of neutrino oscillations renormalizable Supersymmetric SO(IO) GUTs[ll 
[21 m m [5l El [7] have increasingly come to occupy the middle ground of the Grand Unification 
landscape due to their natural accommodation of seesaw mechanisms for neutrino masses 
and their ability to justify R-parity as a part of the gauge symmetry [5] while preserving it 
unbroken to low energies|9l |3] when the most economical and viable Higgs set i.e 126 © 
126 © 210 is used. They thus also predict a stable LSP which is welcome as cosmological 
dark matter. The ability of these theories to accommodate all the known (and probably also 
all the as yet unmeasured) fermion mass-mixing data has also been demonstrated provided 
the full complement of Fermion mass (FM) Higgs allowed by group theory i.e 10, 120, 126 
are useddH [131 [H [IS [16], because the model without the 120 fails^ [IH [12] : thus 
providing a welcome indication that it has matured sufficiently as a scientific hypothesis so 
as to be vulnerable to falsification. 

Given these successes it is natural to ask whether -and with what additional constraints 
on its parameters- this so called New Minimal Susy GUT (NMSGUT)[7] is able to account 
for the longstanding discrepancy between the increasingly accurate experimental estimates 
of the SM gauge couplings at Mz and the values predicted by Renormalization Group 
flows in the NMSGUT. We have shown that the effects of threshold corrections due to 
superheavy particles are modest and compatible with the one loop architecture of MSSM 
coupling unification^ [7]. Moreover we found that the bulk of the viable parameter 
space (according to our 10% error allowances ) imphed that the unification scale(and with 
it all other superheavy masses) was raised by about one order of magnitude or more thus 
alleviating the crisis regarding the too short d = 5 proton decay lifetime estimates in Susy 
GUTs. 

However our analysis of the RG constraints followed a somewhat non standard format. 
We calculated the effects of threshold corrections on the values of sin^ 9w{Ms), ac, Mqut 
rather than the standard choices sin^ 9w{Mz), ac, Mqut (given a^i^Mz)) or a^i^Mz), a^, Mqut 
(given sin^ 9w{Mz)). We concluded that allowing for uncertainties ~ 10% in our knowledge 
of the unification parameters restricted, but by no means eliminated, the parameter space 
of this class of Supersymmetric GUTs. However the difference from the accepted format 
has led to some difficulty in communicating our results to workers in the field. Therefore in 
this letter we carry out a precision analysis in terms of the superheavy threshold corrections 
to the prediction of 0:3 (M^), ac, Mqut- 

Our analysis is aimed at restricting the parameter space of the NMSGUT based on the 
requirement that the parameters be such as to yield a lowering of the prediction of 0:3 (M^) 
from the two loop corrected value of 0.130 which is uncomfortably larger than the central 
value of the world average experimental value 0.120 ± .01. We indicate the region of the 
control parameter space of the NMSGUT (as well as its complex parameter cousin) which 
yields threshold corrections of the right sign and magnitude to bring the prediction for 
as{Mz) back to the observed central value. 
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2 Threshold corrections 



The two loop RG flow equations for the three gauge couplings of the standard model and 
MSSM can be integrated[T7l[l8l[T9] to yield predictions for the Grand Unification scale, the 
value of the gauge coupling at unification and the value of one of the gauge couphngs at the 
scale Mz given the values of the other two as inputs. Since it is a^{Mz) which carries the 
largest uncertainty (~ 8%) while a em{M z), sin^ 9 w{Mz) are quite precisely known (better 
than 0.01%, 0.1% respectively) it is usual [19] to choose to predict a^i^Mz). Using 
updated parameter values [20j 

ll7GeV ; Mz = 91.1876 ± m21GeV 
127.918 ±.018 ; 4 = -23122 ± .00015 
172.7 ± 2.9Ge\/ (1) 

we find from the equations of [19] 

as{Mz) - = 0.130 ± 0.001 + 3.1 x IQ-^GeV-^ x [(mr'')^ - {l72.7GeVf] + H^^ (2) 

where A„ = A^^"^ + Af"**^ threshold corrections. 

The effect of the two loop Yukawa coupling corrections Ha^ was estimated [19] to be 
bounded :-0.003 < H^XKhb) < 

Thus if we place credence on the central value of 0.12 for as(M^), the effects of the 
low energy thresholds ( equivalently parameterized by an effective supersymmetry breaking 
scale MsusY^^)-, GUT thresholds (and even NRO corrections) lumped together into A^^ 
should correct the excessive value found. The effect of the Susy thresholds can raise or lower 
the value of as{Mz). For 2bQGeV > Msusy > 20GeV{the actual values of superpartner 
masses will be much higher since Msusy is a composite parameter formed from all the 
superpartner masses [19]) one finds that 0.005 > Af"*^ > —0.003. Thus it appears that 
as{Mz) — '^a^^ could be as high as 0.135 or as low as 0.124 so that superheavy threshold 
corrections in the range —0.004 > A'^^'^ > —0.015 are indicated. 

As we have described at length in earlier papers[6l|ini[IIl[7], the superheavy thresholds 
are controlled by one 'fast' control parameter x and several other 'slow' parameters which 
are generically of order 1 or else do not affect the thresholds significantly. Thus it is possible 
to scan over the complex x plane calculating the changes in the RG flow predictions due 
to the GUT thresholds and obtain an overview of the behaviour over the whole parameter 
space. By imposing stability of the one loop analysis the viable ranges of x and even of 
the slow parameters can be depicted[6l [HI [7]. In earlier papers we had taken as{Ms) 
as input and restricted changes in the output sin^ 9w (Ms) to be no more than 10%. We 
then obtained a relatively wide range of viable parameter values. With sin^ 9 w{Mz) as 
input and the narrow range of desired threshold corrections as described above one gets 
a relatively narrow range of permitted values of the fast parameter x. However it should 
be remembered that the x parameter is determined by solving a cubic equation with one 
parameter ^ which is a ratio ^ = Xm/rjM of the MSGUT parameters, so that even a fixed 
value of X would allow a ball of the actual parameter values. In [7] we had advocated a 
version of the theory with real parameters only (with CP violation arising spontaneously). 



Mh = 
a{Mz)-' = 
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In that case one needs to scan only over a parameter line rather than the complex plane. In 
the next section we briefly recapitulate the structure of NMSGUT to make the RG analysis 
intelligible. 

The mass formulae required to compute the threshold corrections are all given in [6], [7]. 
We do not repeat them here but simply furnish plots illustrating the restriction of the 
parameter space due to the above described demands placed on it . 



3 Essentials of NMSGUT 

The NMSGUT [7] is a renormalizable globally supersymmetric 5*0(10) GUT whose Higgs 
chiral supermultiplets consist of AM(Adjoint Multiplet) type totally antisymmetric tensors 
: 210{^ijki), T26{T>ijkim), '^'^G(J^ijkim){'i, j = 1...10) which break the GUT symmetry to 
the MSSM, together with Fermion mass (FM) Higgs 10 (Hj) and 120{Oijk). The 126 
plays a dual or AM-FM role since it also enables the generation of realistic charged fermion 
and neutrino masses and mixings (via the Type I and/or Type II Seesaw mechanisms); 
three 16-plets ^'a(^ = 1,2,3) contain the matter including the three conjugate neutrinos 
(i'l)- The superpotential (seejH [211 [5l [6l [7] for comprehensive details ) contains the mass 
parameters 

m : 210^; M : 126 ■ T26; Mh : 10^; mo : 120^ (3) 

and trilinear couplings 



A 

k 

c 



210^; r/ : 210 ■ 126 ■ 126; 7 © 7 : 10 ■ 210 ■ (126 © 126) 

10 120 -210; p : 120 ■ 120 ■ 210 



120 -210 126; C : 120 ■ 210 ■ 126 (4) 

In addition one has two symmetric matrices Hab, Jab of Yukawa couplings of the the 
10, 126 Higgs multiplets to the 16.16 matter bilinears and one antisymmetric matrix qab 
for the coupling of the 120 to two 16 s. It was shown [H] [T5] that with only sponta- 
neous CP violation, i.e with all the superpotential parameters real, it is till possible to 
achieve an accurate flt of all the fermion mass data which furthermore evades the diffi- 
culties encountered in accommodation with the high scale structure of the MSGUT[TT] 
provided [T6| [71 [TH [151 [22] one takes the 10, 120 yukawa couplings to be much larger 
than those of the 126 so that Type I neutrino masses are enhanced. 

The GUT scale vevs and therefore the mass spectrum are all expressible in terms of a 
single complex parameter x which is a solution of the cubic equation 

8x^ - ISx^ + 14a; - 3 = -^(1 - xf (5) 

where £ = — . 

Spontaneous CP violation implies that x must lie [7] on one of the two complex solution 
branches x±{^), (^ G (—27.917, 00)). Since A, 77 are already counted as independent x+(^) 
counts for M/m. 

The other parameters besides x (or equivalently ^) have a much weaker effect on the 
spectra and therefore the threshold corrections and can thus be fixed at representative 
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values ~ 1 while scanning the behavior over the x-plane. Once the viable regions are 
identified the limits of permissible variation in these slow parameters can be explored[7]. 
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Figure 1: The allowed region(other than white) for threshold corrections in the range 
—0.015 < a^{Mz) < —.005 and for which the changes in ac-, Mx are acceptable. The black 
region corresponds to 1 < Ax < 2, and dark grey to < < 1, while the thin light grey 
rim is for —1 < Ax < 0. 



4 Correcting 0:3 (M^) 

In our previous work we had taken it as self evident that, given our demonstration that 
corrections to the unification parameters could be quite modest (~ 10%), there would be 
no problem in accommodating any particular desired values of corrections in this range. 
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( A=1.0, 0.95, y= 0.7, y=0.6, S=.7, £'=.85, k=.75, mo=.8, p=.15) 




Figure 2: Plot of Aa^lMz) against ^ on the CP violating solution branch a;+(^) at rep- 
resentative values of the diagonal parameters. Horizontal lines mark out where Aa^lMz) 
is in the desired range —0.015 < Aa^lMz) < —.005 and determine the range of viable 
(-3.8,-2.5)U(2.4,3)). 



given the number of parameters that affected the superheavy mass spectrum. In other 
words fixing the GUT superpotential parameters for the AM Higgs fields on the basis of 
unification requirements alone seems a vain hope since there is likely to remain a significant 
multi parameter ball compatible with current measured values and limits. However we have 
recently shown[23j that in fact proper numerical fits of the fermion data actually determine 
complete candidate NMSGUT parameter sets. Smoking gun discoveries of supersymmetry 
and GUT characteristic processes such as proton decay will, however, still be necessary 
to choose among candidate fits and thus actually fix the NMSGUT parameters. Here we 
merely restrict ourselves to displaying examples of parameter regions where the threshold 
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corrections are capable of lowering a^lMz) to the current central experimental values while 
maintaining the very accurately known aem{Mz),sm'^ 9w{Mz) at the known central values 
and limiting the variation in ac to 10% and 2 > ALoqiqMx > — 1. Note that although a 
negative change of 10% in LogioMx would be disastrous such a problem in fact never arises : 
in all cases the region allowed by the other parameters always features a raised value of Mx- 
As we have already emphasized previously [7], 122] this raise is in fact a welcome mollification 
of the mounting tension regarding too large d = b operator mediated proton decay rates. 
All RG formulae and mass spectra required for our analysis have already appeared [2T1 El E] 
and will not be repeated here. 



( A=1.0, r]= 0.95, T= 0.7, T=0.6, ^=.7, £'=.85, k=.75, mo=.8, p=.15) 
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Figure 3: Plot of Ax = Alogj^gMx against ^ on the CP violating solution branch 
at representative values of the diagonal parameters. The regions with desired Aas{Mz) ■ 
—0.015 < a3{Mz) < —.005 are marked with vertical pairs of lines and have Mx raised 
above the one loop value. 



6 



In Figure 1 we plot the region of the x— plane which can give threshold corrections 
to a3{Mz)oi the desired size while maintaining acceptable values for a^Mx- As already 
noted in [7j the unification scale is significantly raised over almost the entire allowed range. 
This can be seen in Figure 1 where the region with —0.015 < Aa^lMz) < —0.005 is 
stratified according to the predicted values of Mx- 

In[7] we had identified the contour a;+(^)|^ G [—27.917, oo) in the x plane corresponding 
to a complex solution of the cubic equation for x(with positive imaginary part) as the 
appropriate locus on which to examine the behaviour of a theory in which CP violation 
arose spontaneously even though initial values of all parameters were real. In practice a 
plot for ^ G (—10, 10) is suffcient^. In the present context we wish to indicate which parts 
of the contour allow acceptable threshold corrections for 03 (M^) and the values of ac, Mx 
are within 10% of their one loop values. In Figure 2 we plot the superheavy threshold 
correction to a3{Mz) versus the parameter ^ on the branch x+(,^). In Figure 3 we plot 
A logj^o versus ^. 

From these figures it is apparent that there are significant regions of parameter space 
where the threshold corrections are of the right magnitude and sign to correct the calculated 
2-loop value of 0:3 (M^:) to its observed value. In particular we find that the value of Mx is 
always raised above the one loop unification scale. Moreover this change is accompanied [7] 
by an parallel raising of all superheavy masses so that the dimension 5 operator mediated 
proton decay rates are reduced. 



5 Discussion 

In this letter we have illustrated how not only the general NMSGUT based on the Higgs 
system i.e 126 © 126 © 210 © 120 © 10, but also its more constrained but fully realistic 
version with only real parameters, are comfortably able to provide the threshold corrections 
of the right magnitude to remove the discrepancy between RG corrected and experimental 
a^i^Mz). In fact these constraints from the RG analysis must be confronted [231 El] with 
the determination of GUT parameters that arises whenever a fit of the fermion data ex- 
trapolated to Mx is performed. The question of whether otherwise viable fits which will 
determine all parameters and thus threshold corrections will also keep Aa^lMz) in desired 
range is thus amenable to resolution in coming years. 
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